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We show that homogeneous line broadening drastically affects the performance of atomic Faraday 
filters. We use a computerized optimization algorithm to find the best magnetic field and temper¬ 
ature for Faraday filters with a range of cell lengths. The effect of self-broadening is found to be 
particularly important for short vapour cells, and for ‘wing-type’ filters. Experimentally we realize 
a Faraday filter using a micro-fabricated ®^Rb vapour cell. By modelling the filter spectrum using 
the ElecSus program we show that additional homogeneous line broadening due to the background 
buffer-gas pressure must also be included for an accurate fit. 


I. INTRODUCTION 

Devices utilising thermal atomic vapour cells are of 
increasing interest since they offer high precision with a 
compact and relatively simple apparatus. Examples of 
atomic vapour cell devices include magnetometers [1, 2], 
gyroscopes [3, 4], clocks [5, 6], electric field sensors [7], 
microwave detectors [8, 9] and cameras [10-12], quantum 
memories [13-15], optical isolators [16], laser frequency 
references [17] and narrowband optical notch [18, 19] and 
bandpass filters [20, 21]. 

Making these devices more compact, power efficient 
and lighter is currently a burgeoning area of research [22- 
24], since it allows them to become practical consumer 
products. Particularly for devices that require an applied 
magnetic field, compact vapour cells [25-31] offer the ad¬ 
ditional advantage that small permanent magnets can be 
used to create a uniform magnetic field across the vapour 
cell [32], while consuming no power. However, when con¬ 
fining the atomic vapour in small geometries, additional 
effects may need to be taken into account. For example, 
atom-surface interactions become important for atoms in 
hollow-core fibres [33] or nano-metric thin cells [34]. Also, 
cells with a shorter path length require the medium to 
be heated more to increase the atomic number density. 
Not only will this increased heating cause more Doppler 
broadening but the increased number density will mean 
that self-broadening [35, 36] must be taken into account. 
In this article we investigate the effects of these homo¬ 
geneous and inhomogeneous broadening mechanisms on 
the performance of Faraday filters. 

Faraday filters were proposed in 1956 by Ohman [20] 
for astrophysical observations. They were later applied 
to solar observations [37, 38] and used to frequency sta¬ 
bilize dye lasers [39-41]. In the early 1990s the subject of 
Faraday filters was revived [42, 43]. Such filters have re¬ 
ceived increasing attention ever since, owing to their high 
performance in many applications. Faraday filters now 
find use in remote temperature sensing [44] , atmospheric 
lidar [45-48], diode laser frequency stabilisation [49-51], 
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Doppler velocimetry [38, 52, 53], communications [54] 
and quantum key distribution [55] in free space, optical 
limitation [56], filtering Raman light [57], and quantum 
optics experiments [58, 59]. 

The Faraday-filter spectrum is sensitive to many exper¬ 
imental parameters and so a theoretical model is useful 
for designing filters. However, there are only a few arti¬ 
cles describing computer optimization [60, 61]. In this ar¬ 
ticle we use computer optimization to find the best work¬ 
ing conditions for compact Faraday filters. We find ho¬ 
mogeneous broadening is particularly important for Fara¬ 
day filters in ‘wing’ operation [61, 62] and less so for ‘line- 
centre’ operation [60, 63]. The homogeneous broadening 
mechanism of self-broadening is particularly important 
to include since it is unavoidable at high density. Previ¬ 
ous theoretical treatments of Faraday filters [62, 64, 65] 
have not included the effect of self-broadening; we find 
that self-broadening is important for short cell lengths 
and must be included in the model in order to find the 
best working parameters. The structure of the rest of 
the article is as follows: In section H we introduce the 
typical experimental arrangement for Faraday filters and 
qualitatively explain how they work. In section HI we ex¬ 
plain the computer optimization technique used to find 
the best working parameters and show the importance of 
self-broadening for shorter cells. Section IV describes an 
experiment performed to compare with the theoretical 
optimizations. The results show that buffer gas broaden¬ 
ing and isotopic purity strongly effect the filter spectrum. 
Finally we draw our conclusions in section V. 

II. THEORY AND BACKGROUND 

An atomic Faraday filter is formed by surrounding an 
atomic vapour cell with crossed polarizers (see figure 1). 
When an axial magnetic field (B) is applied across the 
cell, the medium becomes circularly birefringent caus¬ 
ing the plane of polarization to rotate as light traverses 
the cell (the Faraday effect [66]), which leads to some 
transmission through the second polarizer. For a dilute 
atomic medium the effect is negligibly small except near 
resonances, and since atomic resonances are extremely 
narrow, this results in a narrowband filter. If the signal 
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FIG. 1. Illustration of the experimental arrangement. A 
micro-fabricated 1x1x1 mm^ ®^Rb vapour cell is placed 
between two axially magnetized ring magnets. This arrange¬ 
ment is then placed between two crossed polarizers, forming 
the filter. The filter is tested by passing a laser beam through 
and onto a photodiode. The filter transmission is defined as 
the intensity of light transmitted through the second polarizer 
(/j;) divided by the initial intensity before the cell (lo). Light 
out of the passband frequency is either scattered in the cell 
or rejected at the second polarizer (ly). 

being detected is unpolarized then half of the light will 
not pass through the first polarizer. This limits the filter 
transmission to 50%. However using a polarizing beam 
splitter allows one to arrange two Faraday filters to allow 
each polarization component through with little loss [46] . 

In a similar way, if the magnetic field is perpendic¬ 
ular to the light propagation direction, one can also 
make a ‘Voigt filter’ [67] which exploits the Voigt ef¬ 
fect [68]. However, in this paper we will only con¬ 
sider Faraday filters. We have chosen to consider the 
D 2 (n^Si /2 —s- n^P 3 / 2 ) lines of potassium and rubidium 
where n = 4 or 5 respectively. 

For a given cell length the parameters that affect the 
Faraday hlter transmission spectra are the applied held 
(B) and cell temperature (T). The effect of T is predom¬ 
inantly to change the atomic number density [69] and 
secondly Doppler width, while B causes the circular bire¬ 
fringence and dichroism. In general the hlter spectrum is 
a complicated function of these two parameters, due to 
the large number of non-degenerate Zeeman shifted tran¬ 
sitions, each with different transition strengths in which 
their lineshape prohles partially overlap. However, it is 
possible to accurately compute the hlter prohle with a 
computer program [60, 62, 70]. 

We use the ElecSus program to calculate the hlter spec¬ 
trum. The full description of how the program works can 
be found in ref. [70]; here we summarize the key points. 
An atomic Hamiltonian is built up from contributions 
from hyperhne and magnetic interactions. The eigenval¬ 
ues allow the transition frequencies to be calculated while 
the eigenstates can be used to calculate their strengths. 
The electric susceptibility is then calculated by adding 
the appropriate (complex) line-shape at each transition 


frequency, scaled by its strength. The imaginary part 
of these line-shapes have a Voigt prohle [71], which is a 
convolution between inhomogeneous broadening (Gaus¬ 
sian prohle from Doppler broadening) and homogeneous 
broadening (Lorentzian prohle). Typically, the full-width 
half maximum of the Lorentzian has contributions from 
natural broadening (Fq) and self-broadening (Fseif) and 
buffer gas pressures (Fbuf)- The real part of the elec¬ 
tric susceptibility can be used to calculate dispersion, 
whilst the imaginary part can be used to calculate ex¬ 
tinction [72]. This allows the calculation of a variety of 
experimental spectra, of which the Faraday hlter spec¬ 
trum is one. The result is given as a function of global 
detuning, A, which is dehned as A = w — wq, were w 
is the angular frequency of the laser light and wq is the 
global line-centre angular frequency. 

III. OPTIMIZATION 
A. The simple approach 

The optical signal in a vapour cell device comes from 
the interaction of the light with all the atoms in the beam 
path. This means that for compact vapour cells with 
shorter path lengths, the atomic number density must 
increase to compensate for the loss of signal. For ex¬ 
ample the Faraday filter spectrum can be thought of as 
some function of the product crJVL, where A/" is the num¬ 
ber density, L is the length of the medium and a is the 
microscopic atomic cross-section (describing the effect of 
extinction and dispersion due to a single atom). Assum¬ 
ing a remains constant, we can achieve the same filter 
when reducing L by increasing A/" by the same factor. 
Therefore, once good parameters of B and T are found 
for a particular cell length, we can find the new appro¬ 
priate parameters by changing the temperature such that 
AfL remains constant. 

However, this argument will break down at some point 
since a is not generally constant. By increasing the 
cell temperature we also change the amount of Doppler 
broadening. Also, at high densities, interactions between 
atoms cause self-broadening, which can be modelled as 
Fseif = PAf, where /3 is the self-broadening parame¬ 
ter [36] . Both the Doppler and self-broadening will affect 
a. To find where these effects become important we need 
to compare it with a computer optimization technique, 
which can find the best parameters at each cell length. 

B. Computerized optimization procedure 

Efficiently finding the optimal experimental conditions 
for a Faraday filter requires three tools. First a computer 
program is needed which can calculate the spectrum with 
the experimental conditions as parameters. Secondly, a 
definition of a figure of merit (or conversely a ‘cost func¬ 
tion’ [73]) is then needed to numerically quantify which 
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filter spectra are more desirable. Finally, this figure of 
merit is then maximised (or the cost function is min¬ 
imised) by varying the parameters according to some al¬ 
gorithm. 

We used a global minimization technique [74] 
which includes the random-restart hill climbing meta¬ 
algorithm [73] in conjunction with the downhill simplex 
method [75] to find the values of B and T which maxi¬ 
mized our figures of merit. This routine was used in con¬ 
junction with the ElecSus program [70] which calculated 
the filter spectra. ElecSus was used because it includes 
the effect of self-broadening, which is essential for this 
study, and also because it evaluates the filter spectrum 
quickly (< Is) which makes this kind of optimization 
practical, since the filter spectra need to be evaluated a 
few thousand times. 



L (mm) 


C. Figure-of-merit choices 


The signal-to-noise ratio of a narrowband signal in 
broadband noise is greatly improved by using a band¬ 
pass filter. For the case of white noise, the noise power 
is directly proportional to the bandwidth of a top-hat fil¬ 
ter. For a more general filter profile, the equivalent-noise 
bandwidth (ENBW) is a quantity which is inversely pro¬ 
portional to the signal to noise ratio, and is defined as 


ENBW = 


dx iys) 


( 1 ) 


where Ix is the light intensity after the filter, v is the 
optical frequency and Us is the signal frequency. If there 
is freedom in the exact position of the signal frequency 
we can set it to the frequency which gives the maximum 
transmission {Ix{vs) = /max)- 

Although minimising the ENBW is desirable, this usu¬ 
ally comes with a reduction in transmission [60]. Using 
the following figure of merit. 


FOM = 


Ix{v)<lv 


ENBW 


( 2 ) 


) —^ir 


we can maintain a reasonably large transmission [60], 
while minimizing the ENBW. When optimising using 
this figure of merit we often find a wing-type filter spec¬ 
trum [61]. In order to compare with line-centre filters we 
also use the following figure of merit. 


FOM' 


Ix{iy)diy 


/Stt 


( 3 ) 


where we set Vg to be the line-centre frequency. 

To calculate these figure-of-merit values we simulate 
filter spectra with a range of 60 GHz around the atomic 
weighted line-centre with a 10 MHz grid spacing. The 
integration is performed by a simple rectangle method. 


FIG. 2. The figures of merit of filter spectra found by op¬ 
timization or extrapolation. The hollow (olive) circles show 
the figure of merit found by taking the optimal magnetic field 
and temperatnre of the 100 mm length cell and changing the 
temperature such that J^L = const. The solid (purple) dots 
show the figure of merit maximized by changing the magnetic 
field and temperature for each cell length. The main panel 
shows the results of a wing-type filter using an isotopically 
pure ®^Rb vapour, the inset shows a line-centre filter with a 
potassium vapour at natural abundance. Both are modelled 
for the D 2 line of the respective element. 


The limitation to the accuracy of calculated the figure- 
of-merit values comes from the grid spacing; a finer grid 
spacing of 1 MHz only improves the accuracy by 0.2% at 
best. 


D. Results for wing and line-centre filters 


The figure of merit of equation (2) was maximized 
while simulating an isotopically pure ®^Rb vapour with 
L = 100 mm, finding the optimal values of B and T to 
be 67.3 G and 60.9 °G respectively. We then used the 
simple approach (section HI A) to find the new values of 
the vapour cell temperature for a range of shorter cell 
lengths, and then evaluated the figure-of-merit values. 
In addition the figure-of-merit values were re-optimized 
(section HIB) for each cell length to see if further im¬ 
provement could be found. Figure 2 shows the compar¬ 
ison of the two methods. We can see that the figure of 
merit changes with cell length, as is expected, since line 
broadening means that the filter spectra cannot be made 
identical for different cell lengths. We can also see that 
moving to shorter cells has a deleterious effect, but can 
be somewhat mitigated by re-optimization at each cell 
length. 

The inset of figure 2 shows the result of a similar anal¬ 
ysis for a potassium vapour at natural abundance [79], 
this time using the figure of merit of equation (3) to pro- 
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FIG. 3. Atomic number density after computer optimisation 
(Afopt) multiplied by cell length (L), as a function of L. The 
optimisation involves changing cell magnetic field and tem¬ 
perature of the cell in order to maximise the figure of merit 
at each cell length. The dark grey (purple) dots show the 
results when self-broadening is included in the model for the 
filter spectrum, while the light gray (blue) circles show the 
result without self-broadening. The main panel shows results 
for an isotopically pure ®^Rb vapour while the inset gives the 
results for potassium at natural abundance. 



FIG. 5. Gomputer optimized Faraday filter spectra as a func¬ 
tion of linear detuning. The optimal parameters were found 
to be B = 67.3 G and T = 60.9 °C for the 100 mm long ®^Rb 
vapour, B — 85.8 G and T = 127.8 °G for the 1 mm long 
®^Rb vapour, B = 801 G and T = 70.2 °G for the 100 mm 
long K vapour, and B = 864 G and T — 136.1 °C for the 1 
mm long K vapour. The ENBW is 2.0 and 2.2 GHz for the 
®^Rb vapour at 100 and 1 mm length receptively, whereas for 
the K vapour the ENBW is 2.4 and 2.6 GHz at 100 and 1 mm 
length receptively. 



FIG. 4. Filter transmission solid black curve) and cell 

transmission ((B -|- Iy)/Io, dashed blue curve) as a function 
of linear detuning (A/ 27 r), zoomed around the region of peak 
transmission. The left panel models a ®^Rb vapour on the D 2 
line, while the right panel models the K D 2 line; both of length 
1 mm. The cell parameters were set to B = 85.8 G and T — 
127.8 °G {Af = 3.2 x 10“ cm"®) for ®'^Rb, and B = 864 G and 
T = 136.1 °C {Af = 6.0 X 10®^ cm"®) for K. The uppermost 
lines were calculated with a Lorentzian width given by natural 
broadening only (~ 6 MHz) while the middle and lower lines 
have a further 50 and 100 MHz of Lorentzian width. The 
global line-centres occur at 384.23042812 THz [76, 77] for the 
Rb D 2 line and 391.01617854 THz [78] for the K D 2 line. 


duce a line-centre profile filter. The main difference in 
the results is that the figure of merit is less affected by 
decreasing cell length than the wing-type filter. 

The reason for the difference between wing-type and 
line-centre filters can be elucidated by plotting the AfL 
product as a function of L after computerized optimiza¬ 
tion at each cell length, as shown in hgure 3. By repeat¬ 
ing the optimization with the effect of self-broadening 
‘turned off’, we can see that the ®^Rb wing-type filter 
is affected far more by self-broadening than the K line- 
centre hlter. One can understand this difference in the 
behaviour of the two types of Liters by inspection of the 
spectra (see figure 4). Increases in Lorentzian broaden¬ 
ing cause a decrease in transmission through the vapour 
cell at the filter frequency. This happens far more for 
the wing-type than line-centre filters. Changes in trans¬ 
mission on the wing of an absorption resonance due to 
Lorentzian broadening is due to the fact that Gaussian 
broadening decreases much faster than Lorentzian broad¬ 
ening with detuning from resonance [80]. A higher op¬ 
tical depth transition feature will show this effect more 
strongly. This is one of the differences between wing 
and line-centre type filters. Wing-type filters rely on the 
sharp decrease in transmission caused by the atomic reso¬ 
nances to create narrow hlter transparencies. This means 
that the circular dichroism cannot be too large since both 
polarizations need be scattered in the cell to sharply re¬ 
duce the filter transmission to zero. However, a small 
amount of dichroism means that there is a small relative 
birefringence, which means that a high number density is 
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FIG. 6. Experimental and theoretical Faraday-filter spectra on the rubidium D 2 line as a function of linear detuning (A/27r) 
from the weighted line-centre (384.23042812 THz [76, 77]). A 1 mm length vapour cell was used with an isotopic ratio of 99% 
®^Rb to 1% ®®Rb. The solid black line in panel (a) shows the experimental filter spectrum and the dashed (red) line shows the 
fit to theory that includes the natural, self, and buffer gas induced (Fbuf) Lorentzian broadening effects. Below panel (a) the 
residuals, R, (the difference between experiment and theory) are plotted. There is an RMS deviation between experiment and 
theory of 0.6%. The inset of panel (a) shows the effect of Fbuf on transmission (solid purple line) and ENBW (dashed blue line) 
of theoretical filter spectra. The vertical dashed line marks the amount of buffer gas broadening seen in the experiment. Panel 
(b) shows a zoomed in region around the peak at 3.1 GHz, including theoretical curves with natural homogeneous broadening 
only (dashed blue) and with natural and self-broadening (solid blue). 


required to create the large absolute birefringence neces¬ 
sary for the rotation of 7r/2. Conversely, the line-centre 
hlter works by having a large circular dichroism, such 
that the transitions which absorb each polarisation of 
light are almost completely separated. We can see this 
in figure 4 where there the cell transmission is optically 
thick for just one circular polarization on either side of 
the transparency (causing fa 50% transmission of linearly 
polarized light through the cell and « 25% transmis¬ 
sion though the filter). This large dichroism comes with 
a large relative birefringence, meaning that the number 
density can be lower for a line-centre filter. 


Line broadening clearly has a deleterious effect, how¬ 
ever, good filter spectra for shorter vapour cells can be 
found so long as we change both the B and T to re¬ 
optimize the filter. This is shown in hgure 5 where it is 
evident that the optimal filters achieved for a 1 mm cell 
length closely match that at 100 mm length. 


IV. EXPERIMENT 

To compare theory with experiment for a compact cell, 
we used a micro-fabricated 1x1x1 mm^ isotopically en¬ 
riched ®^Rb cell [27] . The isotopic abundance of ®^Rb was 
found by transmission spectroscopy to be (1.00 ±0.02)%, 
in a similar way to that shown in ref. [32]. This iso¬ 
topic impurity affects the hlter spectra, therefore the hl¬ 
ter parameters were optimized taking this into account. 
We found the optimal parameters to be R = 72.0G and 
T = 137.5 °C, which gave a transmission peak at a de¬ 
tuning of 3.1 GHz. 

The experimental Faraday hlter arrangement is illns- 
trated in hgure 1 . The cell was placed in an oven to heat 
the cell near the optimal temperature, while the applied 
axial magnetic held was produced using a pair of perma¬ 
nent ring magnets. The held inhomogeneity across the 
cell was less than 1%. Two crossed Gian-Taylor polar¬ 
izers were placed around the cell to form the hlter. A 
weak-probe [81, 82] beam from an external cavity diode 
laser was focnssed using a lens (not shown in hgure 1) 
with a 30 cm focal length, and was sent throngh the hl¬ 
ter such that the focus was approximately at the location 
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TABLE I. Maximum transmission (Tmax), equivalent-noise 
bandwidth (ENBW) and their ratio (FOM) for a 1 mm long 
isotopically enriched vapour cell. The magnetic field and tem¬ 
perature were 73 G and 138.5°C respectively. The first row 
represents the fit to the experiment shown in figure 6, while 
subsequent rows give the values after certain physical effects 
were removed (theoretically). 


Spectrum 

Umax 

ENBW (GHz) FOM (GHz"!) 

Fit to Experiment 

0.55 

3.0 

0.18 

No buffer gas 

0.77 

2.6 

0.29 

No self-broadening 
or buffer gas 

0.83 

2.6 

0.31 


of the cell. After the filter, the beam was focussed using 
a 5 cm focal length lens onto an amplified photodiode. 
The laser frequency was scanned across the Rb D 2 tran¬ 
sition, and was calibrated using the technique described 
in ref. [83]. 

Panel (a) of Figure 6 shows the experimental filter 
spectrum plotted with a fit to theory using ElecSus [70] . 
The fit parameters were found to be R = 73 G and 
T = 138.5 °C. The first thing to note is that, due to 
the 1% ®®Rb impurity, the peak transmission occurs at 
= 3.1GHz rather than near -1.3 GHz if the cell 
were isotopically pure (see Figure 5). Also, a further 42 
MHz of Lorentzian broadening was added in addition to 
Fq and Fseif, due to the presence of a small quantity of 
background buffer gas in the vapour cell. This value was 
previously measured by transmission spectroscopy to be 
rbuf/27r = (42±l)MHz. Panel (b) ofFigure 6 shows the 
filter spectrum zoomed into the main peak. In addition 
to the experimental and theory fit is the filter spectrum 
for the optimization that did not include the buffer gas 
broadening. We can see that the additional broaden¬ 
ing drastically affects the filter transmission. Also by re¬ 
moving the effect of self-broadening from the theory, we 
again see a larger transmission. Table I quantitatively 
compares the transmission, ENBW and EOM values for 
the curves shown in figure 6. The inset of Panel (a) 
shows the filter transmission at a detuning of 3.1 GHz 
and the ENBW as a function of Fbuf- The transmission 


decreases while the ENBW increases, showing that the 
performance (as measured by the ratio transmission to 
ENBW) of this kind of Faraday filter deteriorates quickly 
with increasing buffer gas pressures. 

The amount of broadening due to buffer gas pressure 
that we observe, typically corresponds to approximately 
1-2 Torr of buffer gas [84, 85]. The fact that this small 
pressure affects the filter spectra by a large amount shows 
that the wing-type Faraday filter spectra are very sensi¬ 
tive to buffer gas pressure. It has previously been shown 
that non-linear Faraday rotation can be a sensitive probe 
of buffer gas pressure [86] , being non-invasive and using 
a simple apparatus. Our results show that it may be pos¬ 
sible to use the linear Faraday effect instead, for which 
it is easier to model the effect of buffer pressure. How¬ 
ever, it is not yet clear if this is more sensitive than using 
transmission spectroscopy [87] . 

V. CONCLUSIONS 

We have described an efficient computerized method to 
optimize the cell magnetic field and temperature for short 
cell length Faraday filters. From theoretical spectra we 
see that wing-type filters in particular are deleteriously 
affected by homogeneous broadening, while line-centre 
filters are less affected. We perform an experiment to 
realise a wing-type filter using a micro-fabricated 1 mm 
length ®^Rb vapour cell, and find excellent agreement 
with theory. While buffer gasses can enhance some sig¬ 
nals using vapour cells [88], they should be kept to a 
minimum in order to achieve the narrowest Faraday fil¬ 
ters with the highest transmission. 
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